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Abstract 



We discuss leptogenesis in a realistic supersymmetric model of inflation with a low re- 
heat temperature 1-10 GeV. The lepton asymmetry is generated by a decaying right 
handed sneutrino, which is produced after inflation during preheating. The infla- 
tionary model is based on a simple variant of the Next-to-Minimal Supersymmetric 
Standard model (NMSSM) which solves the \x problem, called 0NMSSM, where the 
additional singlet <ft plays the role of the inflaton in hybrid (or inverted hybrid) type 
models. The model is invariant under an approximate Peccei-Quinn symmetry which 
also solves the strong CP problem, and leads to an invisible axion with interesting 
cosmological consequences. We show how the baryon number of the universe and the 
nature of cold dark matter are determined by the same parameters controlling the 
strong CP problem, the /i problem and the neutrino masses and mixing angles. 
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1 Introduction 



Leptogenesis is an attractive mechanism which has been proposed to generate the 
observed baryon asymmetry of the Universe (BAU) fl], The mechanism involves 
the out-of-equilibrium decayQ of a heavy right handed neutrino Nr (or sneutrino Nr 
|| ||). The net lepton number L produced in the decay is then reprocessed into baryon 
number B by anomalous (B+L) violating sphaleron interactions, which otherwise 
conserve (B-L) |J . The same physics that allows the right handed neutrinos to decay 
into light leptons is also responsible for a see-saw neutrino mass matrix. Combining 
the see-saw mechanism with the latest experimental data on neutrino masses J7[] seems 
to favour a scale for the right handed neutrino mass Mr in the range 10 7 - 10 14 GeV. 

What mainly distinguishes the different scenarios of leptogenesis that can be found 
in the literature is the production mechanism for the heavy Nr (or Nr). The lepton 
(baryon) asymmetry has to be produced at some early stage in the cosmological 
evolution of the Universe, at some point after inflation ends and before the time 
of nucleosynthesis. In a sense, leptogenesis/baryogenesis is closely related to the 
inflationary dynamics and the post-inflationary reheating era. The right handed 
sneutrino could be itself the inflaton |§, with L generated during the reheating period. 
If this is not the case, and the reheating temperature {Trh) after inflation is larger 
than Mr, the heavy (s)neutrinos can be thermally produced after reheating, and the 
final lepton asymmetry will depend on the out-of-equilibrium conditions at the time 
they decay [[I], §, |j. However, in supersymmetric models such a large T RH may be in 
conflict with the standard bound Trh < 10 9 GeV in order to avoid an overabundance 
of gravitinos ||. On the other hand, if Trh < Mr the right handed (s)neutrinos would 
have to be produced^ by the out-of-equilibrium inflaton decay, either in perturbative 
decays [|12| or by parametric resonance |13| during preheating. In this way, the out- 



of-equilibrium condition for baryogenesis/leptogenesis is automatically satisfied. The 
other two requirements, baryon number violation and C and CP violation, will be 
provided by the sphaleron interactions and complex phases in the neutrino Yukawa 
couplings respectively. 



In this paper we shall extend the model for inflation proposed in Ref. |14| in order 
to include neutrino masses and implement leptogenesis. The inflation model is based 
on the next-to-minimal supersymmetric model and provides an intermediate scale 
solution to the n problem, and the strong CP problem via the Peccei-Quinn mech- 
anism [TJ[]. The supergravity version of the model [jn| solves the n problem via the 



implementation of a no-scale mechanism, provides F-term inflation from the moduli 



1 Models of leptogenesis based instead in the cosmological evolution of flat directions which carry 
lepton number can be found for example in Ref. ^, ^ . 

2 Strictly speaking, the condition Trh < M R does not rule out thermal production of heavy neu- 
trinos. In most of inflationary models reheating is not instantaneous, and the maximum temperature 
rpmax reacne d is usually much larger than T RH (To), and we could have T RH < M R < T max . Like 
in GUT baryogenesis, there could be models where the decay of the inflaton into heavy neutrinos 
may be suppressed or forbidden, but still they could be thermally produced during the long period 
of reheating [O . 
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fields which are stabilised before and after inflation and there is no moduli problem 
or gravitino problem. The inflationary model is of the hybrid type, characterised by a 
not too large scale for the vacuum energy, V(0) 1 ^ 4 ~ 10 8 GeV and a very low reheat- 
ing temperature 0(1 GeV). It is interesting to study leptogenesis within models with 
such a low reheat temperature since in such models thermal production of Majorana 
neutrinos during or after reheating is impossible, so the production mechanism will 
rely on preheating the fields which occurs during the oscillatory period following the 
end of inflation. This model is particularly interesting since the oscillating inflaton 
fields at the end of inflation do not couple directly to the sneutrinos, but only indi- 
rectly via a coupling to the Higgs doublets. Thus the production of sneutrinos during 
preheating is linked also to the production of Higgs scalars, and since the Higgs scalars 
decay into both radiation and neutralinos it becomes possible to relate the relic den- 
sity of the lightest neutralinos to the baryon number of the universe. Relativistic 
axions are also produced during reheating but these are red-shifted away, although 
later on non-relativistic axions are additionally produced by the usual misalignment 
mechanism and these will contribute to cold dark matter. 

The main advantage of studying a realistic supersymmetric particle physics model 
of inflation is that questions such as the nature of cold dark matter and baryogenesis 
via leptogenesis are related, and determined by the same parameters which control 
the particle physics questions of the fi problem, the strong CP problem, and neutrino 
masses although as we shall see there are many uncertainties at present and many of 
our estimates will have errors of one or two orders of magnitude. 

The layout of the rest of the paper is as follows. In section 2 we summarise the 
main properties of the model, and introduce the right handed neutrinos with the 
usual superpotential suitable for neutrino physics and leptogenesis. An estimation 
of the lepton asymmetry produced is given in section 3. Because the right handed 
neutrinos will decay in less than a Hubble time, much before the inflaton has time to 
decay, we need to check how much of the asymmetry survive the reheating era. This 
is done in section 4. In section 5 we present our conclusions. 



2 A model for inflation and leptogenesis. 



The model of inflation we have proposed is based on the superpotential [O]: 

W = XNHiH 2 + K<pN 2 , (1) 

where Hi, H 2 are the Higgs doublets and 0, N are gauge singlets. The superpotential 
is invariant under a U(1)pq Peccei-Quinn symmetry, which is broken during and after 
inflation by the vevs of <fi and N. The vacuum energy V(0), needed during inflation, 
originates from an F-term of the effective sugra theory JRJ. The Higgs doublets play 
no role during inflation, is the inflaton, and N is the second singlet needed to end 
hybrid inflation. Imposing the slow-rolling and COBE constraints for inflation gives 
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the order of magnitude results: 

k ~ 1(T 10 , (0) ~ (N) ~ 10 13 GeV, V(0) 1/4 ~ 10 8 GeV H(0) ~ O(MeV) , (2) 

with iy(0) the Hubble parameter during inflation. In order to have an effective \x term 
in the Higgs sector of the correct order of magnitude, we require A ~ k. The smallness 
of the couplings is accounted for by the use of higher dimensional operators, so the 
superpotential in Eq.|l] should be viewed as an effective superpotential which originates 
from some intermediate scale solution to the \x problem, as discussed elsewhere [14|. 



Due to the smallness of the (effective) couplings, the fields <j> and iV are very long- 
lived, decaying mainly into axions with a decay rate ~ 10~ 17 GeV. Considering 
only the standard perturbative reheating period following inflation, this will give rise 
to a reheating temperature of order a few GeVs, much below the electroweak scale. 

Right-handed majorana neutrinos are introduced in the model with the usual 
superpotential: 

W = M Ri N Ri N Ri + X LRij L i H 2 N Rj , (3) 

written in the eigenstate basis for the N Ri . The right handed neutrino masses and 
Yukawa couplings Xlr has to be such that they reproduced the observed properties 
of the light neutrino spectrum. The recent data from Super-Kamiokande supports 
*-* v T oscillations with sin 2 2# 2 3 > 0.88 and a mass squared splitting Am| 3 ~ 
(1.5 — 5) x 10~ 3 eV 2 . It also favours the large mixing angle solution for solar neutrino 
mixing, with sin 2 29 12 ~ 0.75 and Am 2 2 ~ 2.5 x 10~ 5 eV 2 , although others solutions 
are not excluded. For numerical estimations, and as a working example, we will use 
the results of a recent estimate of all quark and lepton masses and mixing angles based 
on a string- inspired Pati-Salam model flEfl , although we shall only be concerned with 



the leptonic part of this model. Due to the gauged SU(2) R symmetry the model 
predicts three right-handed neutrinos, and the heaviest one is the one associated with 
the third family with a mass of 10 14 GeV. Although this is the heaviest it nevertheless 
plays the dominant role in generating the atmospheric masses and mixing angle, 
due to the Yukawa structure of the model, leading to an automatic neutrino mass 
hierarchy, according to the single right handed neutrino dominance mechanism [17], 



and bi-maximal mixing. In terms of the heaviest right handed scale of the order 
Mr 3 ~ 10 14 GeV, and the Wolfenstein expansion parameter A ~ 0.22, the right 
handed neutrino masses and couplings are given by: 

Mr, ~ 10 14 (A 9 , A 5 , l) GeV ~ (lO 8 , 10 10 , 10 14 ) GeV, (4) 
/A 8 A 4 A 2 \ 

Xlr ~ A 7 A 4 A , (5) 



Note that the lightest right-handed neutrino of mass 10 8 GeV is significantly lighter 
than 10 14 GeV but does not give the dominant contribution to physical neutrino 
masses due to its suppressed Yukawa couplings, although it is light enough to be 
produced through preheating. 
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Given that in the inflationary model T RH < 1/(0) 1/4 < M R . } the Nr. (N R J fields 
cannot be produced thermally at any stage, nor in the perturbative decay of and 
N fields. Therefore, we depend upon preheating for that. 

Due to the lack in the superpotential of a direct Yukawa coupling between the 
singlets superfields 0, iV and N R , preheating of the right handed neutrinos does not 
look possible. The situation is different for the scalar components, because of the 
coupling X LR with the doublet Higgs H 2 and the leptons Lj. This will induce a term 
in the scalar potential of the form: 

V = ...+ \\ LR \%\H 2 \ 2 N R N Rk + .... (6) 

Large oscillations in the Higgs fields will be induced through their coupling to the 
singlets in Eq.|I|, and in turn we expect this to trigger the preheating of the sneutrinos. 
The lepton asymmetry will be generated by the decay of the sneutrinos instead of 
that of the neutrinos. We notice also that the maximum possible T that we can 
reach during reheating is going to be smaller than the M Ri masses (at most T max ~ 
10 8 GeV ~ M Rl ), so once produced, the sneutrinos will remain out of equilibrium 
and will decay faster than say the inflaton. 



3 Lepton/Baryon asymmetry 

Preheating of right-handed sneutrinos provides the seed we need for leptogenesis, i.e., 
a non-zero number density of the order 



V(0) 

U ^ R i ~ ~M R ~ ' (7) 



where q parametrize the fraction of the total vacuum energy which is transfered to the 
sneutrinos during preheating. Given the hierarchy in masses, it is not unreasonable 
to assume C3 -C c 2 , c\, whilst we will take c 2 ~ c\ ~ c ~ 0(1 — 0.1). CP violation 
in the decay of N Ri comes from the interference between the tree-level and one-loop 
amplitudes [H| [TDfl . The CP asymmetries given by the interference with the 

one-loop vertex amplitude are [0, |J: 

T(N Ri -> f + gg) - T(N Ri -> g + Hj) 
T(N Ri ^!+H 2 ) + T(N Ri ^d + Hl) 

= gMt \ : E(^[(4^) 1j ] 2 )/(M^./M^), (8) 



where 



l-(l + a:)log(i±^) 



(9) 



The interference with the absorptive part of the one-loop self-energy also gives a 
contribution to the asymmetry, which in general is the same order as those given 
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above, unless the (s)neutrinos were almost degenerate in which case could it be much 
larger @, ITT 



As an example in order to estimate the values of e,, we will consider the model 
given in Eqs. |] and |5|. Assuming maximal CP violation ( Im[- • -] 2 ~ | ■ ■ • | 2 ), we 
can see that the asymmetries will be dominated by the larger couplings to the third 
generation of leptons, with: 

\^lr\1 3 M Ri ^ 9 m-7 1n -8 
ei ~ — —j— ~ 10 - 10 , (10) 

07T Mr 3 8tt 

[Azfilig m R2 a 4 4 5 

e 3 ~ \hl^kM^ „ ^! „ 10-10 _ 10 -u . (12) 
8tt Mr 3 8tt v ; 

We remark again that these are only order of magnitudes estimations, with large 
uncertainties in their values. Besides, the values of e« are model dependent. Other 
texture models with values of the Yukawas conistent with the experimental data on 
neutrinos, and similar hierarchy among the right handed neutrino masses, could give 
rise to a larger asymmetry such as e\ ~ 10~ 6 [2H . 



The decay of the sneutrinos occurs fast enough to neglect any effect due to the 
expansion of the Universe, and the lepton asymmetry is then given by: 

n B -L - ^n Ri ~ e ^TZ^ • ( 13 ) 

The decay of the heaviest right handed sneutrino (if produced) will give rise to a 
negligible lepton asymmetry, whilst that generated in the decay of N Rl and Nr 2 
are comparable because ei /Mr 1 ~ ^/Mr 2 ~ 1/Mr 3 . This is then converted into 
baryon number by B+L violating sphalerons interactions (which are in equilibrium 
for temperatures in the interval ~ [200, 10 12 ] GeV) [PHI , 



n B 
s 



UB = ~23 UB ~ L ' ^ 
and finally at the time of nucleosynthesis we will have: 

8 (cl/(0))V4 7 (cV(0))V4 

~ ^ lei ~^T ^d A ^33 ' , (15) 

nud. 2.6 M Ri ZSn M R3 

where we have used s = (27r 2 /45)g*T 3 evaluated at the time the leptons are produced, 
i.e., T ~ .3(cy(0)) 1//4 , with the effective number of relativistic degrees of freedom 
g* ~ 100. The factor 7 accounts for the dilution due to possible entropy production 
during reheating. Substituting the values of V(0) and Mr 3 Eqs. (0) and (f|), with 
Xlr 33 — 1, we obtain: 

^ -c^xKTV (16) 

S nucl. 

In order to explain the observed baryon asymmetry %/s ~ 10~ 10 []22j we cannot 
allowed much entropy (radiation) production during the reheating era. In the next, 
we will try to estimate the factor 7 based on simple assumptions. 
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4 Preheating/Reheating 



During the oscillations of the background fields (j) and N, particles can be produced 



by parametric resonance (preheating) |23|, 21. 25] much before the inflaton has time 



to decay perturbatively, being this in general a more efficient mechanism of particle 
production than standard perturbative decay. 

Even if the couplings are very small, the amplitude of the oscillations is large 
enough to "preheat" the modes of the scalar fields and N. Beside, the value of 
Hubble parameter is small in the model, which allows for a large number of oscillations 
in a Hubble time, before they start to feel the effect of the expansion. Because hybrid 
inflation ends in a phase transition, with the effective squared mass of the N field 
changing sign, production of (j) and N quanta is very efficient during the first few 
oscillations of the background fields |26|] . Due to the coupling A between the Higgses 
and the N field, we also expect to preheat the Higgs fields Hi in a similar way. The 
evolution equations for the Higgs quantum fluctuations are indeed analogous to those 
of the singlets, and we can assume similar number densities for both. The sneutrino 
fields Nptt are therefore preheated through the Higgses, with the lightest one more 
likely to be produced, and they will be clearly out of equilibrium. We may also 
preheat axions, fermions, etc.. but with much smaller number densities^ 

Preheating is efficient only in producing very low frequency modes. Nevertheless, 
rescattering effects will allow to excite higher frequency modes and redistribute the 
energy density. Based on the results for the singlets |[26|| , we may estimate that after 
just 3-4 oscillations a fraction of the vacuum energy has been transfered to the singlets 
and the other fields quanta, with more or less equal energy densities. The typical time 
scale for this to happen is given by: 

At preh ^^^~10-*GeV-\ (17) 

where AN osc counts the number of oscillations, and is the mass of the fields (j) and 
N in the global minimum (and therefore the typical frequency of their oscillations): 

M = k(N) = 0(1 TeV) . (18) 

The decay rates of the fields involved, Nr 12 , Hi, and singlets and N, can be 
estimated as: 

r N Rl * l -^^M Rl nO(0M)GeV, (19) 



I \ I 2 



Nr 2 



M R2 « 0{ltf)GeV, (20) 



S7T 

T Hi ~ f-MtfW 0(10) GeV, (21) 
^ ^M^O(10- 17 )GeV, (22) 



5 Due to the smalleness of the couplings, non-thermal production of gravitinos |27fl is not a problem 
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where g is the electroweak coupling constant, and the mass of the Higgses are M#. ~ 
0(100 GeV — ITeV). The sneutrinos N R2 will tend to decay immediately after they 
are produced, r^. <C At pre h, and its decay products quickly thermalize by scattering 

from each other given that: 

An t a sc > H (23) 

where An/ is the number density of the light degrees of freedom Am ~ , a sc oc 
M^ 2 2 , and H is the Hubble parameter of the order O(MeV). No back reaction is 
expected from them, except that part of the vacuum energy is converted into radiation 
with a temperature T ~ 0.3c2 /// V(0) 1//4 . The fields Hi and N Rl do not decay before 
At pre h. Therefore, back reaction effects due not only to the singlets but also to 
the Higgses (and eventually N Rl ) will soon slow down and suppressed the rate of 
production of particles during preheating, in less than a Hubble time if -1 . 

At this point we can consider that the Universe has been reheated up to a temper- 
ature T ~ O(10 8 GeV), but with a non-negligible fraction of the energy still in the 
form of cold oscillations and singlets, Higgses and N Rl . The right handed sneutrinos 
N Rl also decay out-of-equilibrium in a time 8t ~ T^ 1 <C H~\ transferring its energy 

to the thermal bath. On the other hand, the decay rate of the Higgses will now 
be suppressed by the factor Mh/T , rendering it quite inefficient. Therefore, after 
preheating becomes inefficient, back reaction and rescattering effects take place, and 
both Nr 2 and N Rl decay, we are left with the vacuum energy distributed among the 
singlets, Higgses, and radiation. 

Preheating can become very inefficient, but will not neccessarily stop as long as 
and N continue to oscillate (that is, there is some energy density left in these 
fields), and production of Higgses from the singlets might continue in a very narrow 
resonance regime ||23|| , at the same time than standard reheating. On the other hand, 
we do not expect this effect to be enough to further induced the production of the 
much heavier right handed sneutrinos. The preheating/reheating era will end when 
the singlets finally decay, at a time t ~ FT . 

If the masses of the singlets are smaller than the lightest Higgs mass, the fields <fi 
and N decay predominantly into axions, with the branching ratio into others particles 
being much smaller. The axions behave as relativistic particles, but they do not 
thermalize [BOJ, that is, their interaction rate always remains smaller than the Hubble 



expansion parameter. The axion interaction rate Tai is given by: 

Tai = Wa\v\)n R (24) 

where n R is the radiation number density n R ~ T 3 /tc 2 , and a a is the axion cross 
section for scattering off the thermalize radiation. On dimensional grounds, the cross 
section can be written as[| a a = a a J f%, f a being the axion decay constant 10 12 — 
10 13 GeV. 



4 The coupling a a is either due to tree-level interactions, and therefore further suppressed by a 
factor T~ 2 when the temperature is larger than the typical mass scale of the particle exchanged, or 
due to loops effects, which are suppressed by a factor l/(87r 2 ) 2 . Therefore, until the temperature 
drops near 0(1 GeV) we have a a ~ 10~ 4 . 
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Let us define to as the initial time after the initial burst of radiation produced by 
the sneutrinos decay. At this time, the Universe is at a temperature T ~ O(10 8 GeV) 
but still H ~ V r (0) 1/2 /v / 3M P . The ratio Y AI /H is then: 



AI 



H 



a a T 3 V3M P _ , /4 f 10 12 GeV\ 2 , 

^f^VW-^ Wa ' c {-^) <u (25) 



After the Universe becomes radiation dominated at a temperature Trx> we will have 
instead: 



AI 



H 



to 




10 12 GeV K 2 



<10- 5 a a ( X " p v ) T RD <1, (26) 



where we have used H ~ (ir 2 gT/90)T 2 /Mp, and g-r is the effective number of relativis- 
ts degrees of freedom at T. The last inequality in Eq. (fffi) follows because the factor 
a a becomes 0(1) only when T ~ 0(1 GeV). If radiation is being produced while the 
Universe is matter dominated, the number density and Hubble parameter evolve as 
tir oc a~ x l 2 and H oc a~ 3 / 2 , with a the scale factor, so the ratio Yai/H increases 
in time. In a radiation dominated Universe they will both scale as a _1//2 . In either 
case, Eqs. ( p5| ) and ( pE|) ensure that the axions never come into equilibrium. This 
means that we do not expect radiation to be produced from the singlets (inflaton), 
contrary to most models of inflation. Any extra radiation will come eventually from 
the Higgses decay. 

Let us now briefly summarise the above discussion. We have argued that right- 
handed sneutrinos are produced during the initial period of preheating, but they 
decay rapidly into leptons and Higgses. We are now interested in the evolution of 
the Universe from this time until the time of reheating defined as the time at 
which the singlets completely decayQ. In this interval there is an interplay between 
the energy density of the oscillating inflaton fields p^, the number density of Higgses 
n,H, the energy of the axion fields p a xion, and the energy density in radiation pp. We 
shall model this as follows, p^ will be steadily reduced due to continual production of 
axions (through the standard perturbative rate T^) and Higgses (through inefficient 
preheating). The axions behave as relativistic matter, but they stay out of equi- 
librium. Higgses are created and annihilate into radiation with a thermal-averaged 
crossed section (an\v\) and decay rate (Tjj). Therefore, the radiation density pr re- 
ceive contributions from the Higgses but not from the singlets. On the other hand, 
among the decay products of the Higgses we will find also neutralinos, the lightest of 
them being a candidate for cold dark matter. The standard calculation of their relic 
abundance |31] depends mostly on their freeze-out temperature, the temperature at 



which they decouple from the plasma, and it is usually assumed that this happens 
while the Universe is radiation dominated. That is the case when the reheating tem- 
perature is much larger than the typical mass scale of the particle, that is, reheating 
ended much before they freeze-out. However, as shown in Ref. P2j, the situation 



5 In general, the reheating time will coincide with the time at which the Universe becomes radiation 
dominated iflD. However in the present context this may not be the case, see below. 
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changes in an scenario with a low reheating temperature, such that freeze-out takes 
place when the Universe is still matter dominated. In fact it is shown that, for a 
B-ino like lightest neutralino, the cosmological constraints on the B-ino mass and/or 
the right handed slepton mass are relaxed and even disappear once a reheating tem- 
perature below the B-ino mass is allowed. In our type of scenario, moreover, the 
Higgses will be kept for a while out-of-equilibrium, due to preheating, and therefore 
neutralinos will be produced also out-of-equilibrium, initially most likely in a matter 
dominated universe. If they do not re-enter equilibrium before reheating is complete, 
their relic abundance may be different than that obtained in other scenarios. 

Under the above assumptions, the evolution of the energy densities (singlets, ax- 
ions and radiation) and number densities (Higgses and neutralinos) during reheating 
can be described by a simple set of equations |L0], [II], [3^| : 



P4> — — "iHp^ — T^p^ — TprehPcf, , (27) 
Paxion AH Paxion ^<f>aP<t> > (^^) 



TV* 

h H = -3Hn H - (1 - B x )(T H )(n H - n e g) - B x (T H )(n H ^ff) 

n x 

-(a H \v\)(n 2 H -nf)+T preh ^f, (29) 

h x = -3Hn x + B x (V H )(n H -^n^)-(a x \v\)(nl-nf), (30) 

n x 

p R = -4Hp R +(l-B x )(r H )(E H ){n H -n%) 

+2(a H \v\)(E H )(n 2 H - nf) + 2(a x \v\)(E x )(n x - nf) , (31) 

where T pre h models the rate of production of Higgses through "inefficient" preheating, 
(E H ) and (E x ) are the average energy per Higgs and neutralino respectively, with 



(Ej-y/mt + T* ,i = H, X , (32) 

and in the numerical calculations we have set the masses m# = = 10m x = 1 TeV, 
and 

(In) ^ a H m H (l - exp{-m H /T)) , (33) 
(*M)*£{\-exp{-T* (34) 

with ctff = a x = 10~ 3 . The preheating period would set the initial conditions to solve 
these equations, and at t = we will take: 

P<t>^ Pr^ (E H )n H ~ . (35) 

When reheating ends at t = tjui, the temperature of the Universe will be Trh oc 
yT^Mp, where Mp is the reduced Planck mass and 1^ = T$ a + Y pre h the total decay 
rate of the singlets. 
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To start with, the Higgses will soon approach equilibrium due to scattering process, 
with (a H \v\)n e H > H(t ). If T preh <C T^, no appreciable amount of Higgses/radiation 
is further produced and the Higgses will later decay in equilibrium. This means that 
pn will be redshifted as a -4 , faster than the axions which are produced in singlets 
decay, p ax ion oc a -3 / 2 , and when reheating ends we will have an axion-dominated 
Universe. Given that the entropy would be conserved, and that PaxionitRn) — P^rh), 
one can estimated the final ratio of the energy density of axions to radiation as: 



Paxion 



Pr(T rh ) p R (t ) \T RH 



10 8 . (36) 



<1.8. (37) 

nucl. 



The relativistic axions will behave as an extra generation of neutrinos at the time 
of nucleosynthesis, with the number of extra generations constraint by 8N V < 1.8 
f22|| . This translates into a bound for the energy density of axions relative to that in 
radiation: 

^ Paxio'i 
%~PR~ 

Comparing this bound with Eq. (|36"D, it is clear that we need the Higgses out of 
equilibrium and start decaying into radiation at some time t\ much before reheating 
is complete. This will be achieved taking T pre h > T,p a . This also means that allowing 
the standard perturbative decay of the singlets into Higgses does not solve the problem 
of axion dominance, since one expects the respective branching ratios into axions and 
Higgses to be at most of the same order. We need some extra effect apart from 
"perturbative" decay. However, if we take T pre h as a constant parameter, and larger 
than T^a, we are forcing the singlets to decay completely through preheating, which 
in a realistic scenario is unlikely to happen. It is more reasonable to consider this 
parameter as a decreasing function of time: as the energy density of the oscillating 
singlets decreases, the rate of production through preheating will also diminish. We 
will consider a simple ansatz where T pre h follows an exponential law, 

^preh{t) = r p reh(^o) e U > (38) 

with b a constant such that b < r prefe (t ) _1 - We can imagine this like T preh (t) be- 
ing switched on for a while, until the energy density p<p diminishes enough to make 
parametric resonance completely negligible at t > ft -1 , when T pre h(t) is switched off. 

The situation now is as follows: the Higgses start decaying at a time t\ < tun, 
which can be estimated as: 



tr 1 - m) - M '!\ Tpr f° ] - w-^ptM , (39) 

n H (ti) M H M H 



and given that until this time the entropy is constant, and T oc a 1 , we have: 

H » Ti' 3 T%l h ~ (l0 2 °rv„„) 2/3 GeV . (40) 

Higgs/radiation production stops at t 2 — r pre / l (to)~ 1 - Inmediately after the Higgses 
decay. Therefore, entropy is only released between T\ and T2, and in that interval the 
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Figure 1: Evolution of the energy densities of the singlets (cf>, N), axions, Higgses, radiation 
(dashed line) and neutralinos (long-dashed line), when T pre h(to) — W~ 12 GeV, and b — lO 5 r0. We 
have taken B x — 0.5. We have also included the factor 7 - " 1 = Sf/Si (dotted line). 

ratio paxionl Pr is frozen. After t 2 we are left still with the singlets producing more 
axions, so the ratio of axions to radiation again increases until tns when is frozen 
again. 

This can be seen in Fig. II where we have plotted the evolution of the energy 
densities for the case T preh (t ) = IO -12 > and b = 10 5 r^ as an example. The 
time scale is given in units of T^. The singlets start decaying through T preh at t ~ t\, 
but still we have assumed that they initially dominated the energy density of the 
Universe. The Higgses initially are in equilibrium, but due to the contribution from 
T preh they start to decay into radiation and neutralinos. At this point the ratio Sf/Si 
starts to increase as can be seen in the figure. Entropy production stops at t 2 , when 
the singlets energy density is partially depleted due to T pre h. The Universe becomes 
radiation dominated at a temperature T 2 ~ Trd — 10 3 GeV; this is larger than the 
lower reheating temperature we would expect only from T^. Soon after, inefficient 
preheating is swichted off and the Higgses decay. Until T^t ~ 0(1) we still produced 
more axions from the remaining singlets. The final ratio of the axions to radiation is 
Paxion/ Pr ~ 10~ 3 , which is consistent with nucleosynthesis. The neutralinos follow the 
same evolution than the Higgses as far as both are relativistic. Once the Higgses decay, 
they go into equilibrium and when they become non-relativistic the ratio n x /s freezes 
out. At T RH we have {n x / s)\rh — 4 x 10 -11 , which would imply a relic abundance of 
neutralinos of order 1 today; they would dominate the dark matter in the Universe. 
We notice that in this example the neutralinos enter into equilibrium and freeze out 
in a radiation dominated Universe, so the calculation of their relic abundance would 
not differ from the standard one, and the usual bounds would apply. 
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In the above example, the Universe becomes radiation dominated at a temperature 
T 2 < T\. However, from Eq. ( flQ|) we see that had we taken T preh (t ) > 10~ 8 , then 
T\ ~ To, and the Universe would be radiation dominated indeed during the whole 
period of what we call "reheating" , that to say, until the singlets completely disappear. 

With this in mind, we now turn to the calculation of the entropy dilution factor 

HUH!) 3 © 3 ' 

which will depend on the initial value of I^e^o), but not on its time dependence. If 
Fprehito) < 10~ 8 , then the value of T 2 is given by the condition p^fo) — Pnih), and 
we obtain: 

7" 1 « 0.4 (%) 3 « / T ° < 10 8 . (42) 

\Tj p R (t 2 ) JT preh{to) M P 



The smaller the value of T pre h(to), the larger the dilution becomes. On the other 
hand, for large values of T pre h(to), entropy production takes place when the Universe 
is already radiation dominated, the radiation energy density scales as oc a -1 , and the 
time £ 2 is given by the condition p</,(£ 2 ) ~ Pnih) instead. In this case, 7 tends to a 
constant value, given by: 






10 4 x _ ( 4 3) 



Due to entropy production, the final ratio of axions to radiation given in Eq. ( p6[) 
is also diluted, such that: 

Paxion(T RH ) _ /y(£p) / Tp \ ^ T preh (t ) ^^ . 

Pr{Trh) Pn(t ) \TrhJ " b 

The last factor in the above equation is due to the partial depletion in the energy 
density of the singlets around the time £ ~ r pre ft(£o) _1 , with: 

F[x] = x (1 - exp(l/x)) . (45) 



To summarize, in order to avoid axion dominance at the end of the reheating 
period, we have allowed the singlets to decay into Higgses through inefficient pre- 
heating, which will decay into radiation. The effect on the final ratio of axions to 
radiation is twofold: on one hand it reduces the ratio because part of the singlets 
have been converted into Higgses instead of axions; on the other hand, the final ratio 
gets also diluted by a factor 7 due to entropy production. However, the same factor 
7 will dilute the initial lepton assymmetry produced in the s-neutrinos decays, which 
was our main concern. The lower dilution factor is obtained when what we have 
called "reheating" starts directly with a radiation dominated Universe instead of the 
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usual matter-inflaton dominated Universe. Based on simplest assumptions, we have 
obtained an upper bound on 7, 



7 ~l(T 4 x^M, (46) 

which means that the baryon asymmetry at the end of reheating will be at most: 

^ ~ 10- 12±2 , (47) 
s 

a couple of orders of magnitude below the observational data. Given that we have 
only considered a kind of toy model to study reheating, and the uncertainties in it, 
we regard this result as quite promising. In particular, we remark again that the 
value of r pre h is an unknown in the model, controlled by the physics of preheating. 
Moreover, we have taken pR(t ) ~ p^(t ), but it may happen that we could produce 
more radiation than expected before reheating starts. Finally, we have considered 
a particular model for neutrinos masses as a working example, such that the values 
of the Yukawa couplings generate an asymmetry of the order e\ ~ 10 -8 . Other 
texture models with values of the Yukawas consistent with the experimental data on 
neutrinos, and similar hierarchy among the right handed neutrinos, could give rise to 



a larger asymmetry such as €\ ~ 10 p0| , which would give then the correct order 
of magnitude for the baryon asymmetry. 

In any case, the less dilution we can have, the better, and this translates into the 
Universe becoming radiation dominated as soon as possible. In general, this means 
that neutralino freeze out will take place in a radiation dominated Universe, and the 
standard bound on their relic abundance will apply. 

Finally, we mention that we could also avoid axion dominance allowing the singlets 
to decay into other light degrees of freedom apart from axions. However, if the singlets 



decay into radiation this will imply a too large dilution factor, with [10 



TsMp 

7~ V , n , -10" 8 , (48) 

MO) 

where the entropy density is s ~ Tq , and the energy density in singlets p^ are taken 
as the initial values at to- 



5 Conclusion 

In this paper we have discussed a realistic supersymmetric model of inflation WA 



which couples the inflaton to the Higgs, and when enlarged to include right-handed 
neutrinos allows leptogenesis as the mechanism to generate the observed baryon asym- 
metry of the Universe. From the particle physics point of view, the model is an exten- 
sion of the Next-to-Minimal Supersymmetric Model and solves the p problem via an 
intermediate scale which generates the vev for the singlets. The interaction between 
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the singlets and matter fields are dictated by an approximate U(l) Peccei-Quinn 
symmetry, providing also a solution for the strong CP problem, and the axions as a 
candidate for dark matter. As a hybrid inflationary model, it has a quite low scale 
for inflation (O(10 8 GeV), and its predicts an spectral index n = 1 consistent with 
the recent Boomerang and Maxima-1 data |33[, and in principle a very low reheating 



temperature of the order of a few GeV, barring the possibility of both GUT and elec- 
troweak baryogenesis. Therefore we have extended the model to include right-handed 
neutrinos, and have appealed to preheating to produce the lightest right-handed sneu- 
trino so that baryogenesis may proceed via leptogenesis. 

Extending the model to include heavy right handed neutrinos is one of the pre- 
ferred solutions not only for the sake of leptogenesis, but in order to generate a light 
neutrino mass spectrum through the see-saw mechanism, given the strong evidence 
from experiments in support of such light masses and mixings. In this paper we 
have used as an example a realistic model of all quark and lepton masses and mixing 
angles |16|] based on single right handed neutrino dominance [|17|]. We have chosen 



this particular example because it predicts not only a hierarchical spectrum for the 
light neutrinos but also for the heavy right handed neutrinos/sneutrino and in par- 
ticular involves a relatively light right handed state which is available for preheating, 
although we emphasise that this spectrum came out of an analysis of neutrino masses 
and mixing angles which was not performed with leptogenesis in mind. The right 
handed neutrinos and sneutrinos are too heavy to be produced thermally at any 
stage in our inflationary model. However, we have argued that at least the lightest 
right handed sneutrino can be produced through parametric resonance during pre- 
heating, due to its couplings to the Higgs doublets. Its CP-violating decay gives rise 
to a lepton asymmetry, later converted by sphalerons into baryon number. Note that 
the right handed sneutrinos do not couple directly to the oscillating inflaton fields, 
but rather indirectly via the Higgs doublets. Thus Higgs scalars are also expected to 
be produced during reheating, and since these decay into radiation and neutralinos 
we may estimate the amount of neutralinos and entropy that is produced during pre- 
heating as shown in Figure |T|. Note that relativistic axions are also produced during 
reheating but these are red-shifted away, although later on non-relativistic axions are 
additionally produced by the usual misalignment mechanism and these will contribute 
to cold dark matter. 

Because the heavy sneutrinos decay long before reheating is completed, the lep- 
ton asymmetry will be subsequently diluted by the entropy produced in the decays of 
singlets and Higgses, before the time of nucleosynthesis. In order to avoid too much 
dilution, we first required the singlets to decay only into axions, which do not ther- 
malize and do not contribute to the radiation energy. The radiation energy density 
has its origin in the out-of-equilibrium decay of Higgses and sneutrinos, which have 
been previously produced during preheating. If no more radiation is produced, by 
the time the singlets completely decay the Universe becomes axion-dominated, vio- 
lating by many orders of the magnitude the bound on the number of extra relativistic 
neutrino-like species at nucleosynthesis. We argue then than "inefficient" preheating 
of the Higgses is required in order to allow their out-of-equilibrium decay and some 
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extra production of radiation. Due to the long life-time of the singlets, it is possible 
that, while they are oscillating, they could "preheat" other fields, at a rate similar 
or even larger than the perturbative decay rate. To illustrate this point we have 
presented a simplified analyses of the "reheating" period, parametrising "inefficient" 
preheating by a rate T pre h, in order to compare it with the perturbative decay rate 
T^, and it is the result of this simplified analysis which is presented in Figure [|. The 
main qualitative conclusion is that in order to avoid axion dominance we would re- 
quire T pre h > T^. Therefore, the main parameter controlling the analyses is the ratio 
of these decay rates (the perturbative one and through inefficient preheating), which 
will also depend to some extent on the ratio of the Higgs and singlets masses. Note 
that the lighter the Higgses are, the easier they will be produced. However, this may 
modify the period of preheating following the end of inflation. The more radiation 
which is produced the less the dilution we will have later (Eq. (|46|) . 

To summarise, we have presented a semi- quantitative scenario for leptogenesis 
in the context of a realistic supersymmetric low scale hybrid inflationary model. A 
novelty of the model is that during the reheating period, the dilution of the lepton- 
baryon asymmetry is not due to entropy produced in the inflaton decays but due to 
Higgs decay. This in turn is controlled by how many Higgses we are able to "preheat" 
from the singlets before reheating is finally completed, and it will help to avoid an 
axion-dominated Universe at the end of reheating. In a model with a hierarchy in the 
masses of the heavy right handed neutrinos, the combination of the small asymmetry 
e with some later dilution could give rise to the correct order of magnitude for the 
final value of ns/s. Note that in a realistic supersymmetric model such as this the 
value of the baryon number is related to the question of the nature and abundance 
of cold dark matter, and that these questions are in turn related to the questions of 
the n problem, the strong CP problem, and neutrino masses. 
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